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Abstract 

(N ! We employ the techniques introduced in the companion papers [l][2|[3] to derive a connection 
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formulation of Lorentzian General Relativity coupled to Dirac fermions in dimensions D+ 1 > 3 
with compact gauge group. The technique that accomplishes that is similar to the one that 
' has been introduced in 3 + 1 dimensions already: First one performs a canonical analysis of 

CO . Lorentzian General Relativity using the time gauge and then introduces an extension of the phase 

ly-^ ' space analogous to the one employed in ffj to obtain a connection theory with SO(D + 1) as the 

internal gauge group subject to additional constraints. The success of this method rests heavily 
on the strong similarity of the Lorentzian and Euclidean Clifford algebras. A quantisation of the 
Hamiltonian constraint is provided. 
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1 Introduction 



In the previous papers of this series [TJ [21 [3] we introduced a new connection formulation of 
vacuum General Relativity with compact gauge group in any spacetime dimension D + 1 > 3. 
In this paper we consider coupling of this theory to matter. Here it will suffice to consider 
fermionic matter because gauge bosons can be coupled in the same way as it has been done 
in 3 + 1 dimensions already [5J: Scanning through the details of [3J one realises that nothing 
depends substantially on D = 3 and since the background independent LQG representation for 
theories of connections as well as its uniqueness generalises to any dimension, we can consider 
the gauge boson sector as treated already. 

Our starting point is the standard canonical treatment of fermions coupled to General Rel- 
ativity: 

To the best of our knowledge, Kibble [5] was the first to consider the canonical formulation of 
fermions coupled to vierbein gravity. The classical coupling of fermions to the new variables [6] 
was provided in [7]. Since then, several papers appeared debating issues arising when including 
fermions. Among others, the role of the Immirzi parameter [8], the appearance of torsion [9|ll0j 
and the correct form of the Hoist modification are ongoing debates. Here we will consider 
the simplest possibility, namely the standard coupling of Dirac fermions to vielbein gravity. 

In 3 + 1 dimensions the quantisation of this theory was carried out for the first time in the 
context of the new variables in [JJ [12] . The new ingredient was the passage to Grafimann valued 
half densities and a representation in terms of holomorphic wave functions of the fermionic 
variables. Technically, in 3 + 1 dimensions one works in the time gauge and with the Ashtekar 
- Barbero connection which can be obtained by an extension of the ADM phase space subject 
to an SO (3) Gaufi constraint. 

In higher dimensions, an Ashtekar - Barbero like connection is not available and therefore 
a new idea is needed in order to arrive at a connection formulation with compact gauge group 
although we are considering Lorentzian gravity. We start from the usual Dirac - Palatini La- 
grangian for Lorentzian General Relativity and introduce the time gauge. This results in a 
formulation in terms of a canonical pair (K l a , Ef), a, b, c, .. = 1, .., D; i,j, k, .. = 1, .., D which 
is subject to an SO(D) Gaufi constraint. We now extend this phase space by a canonical pair 
(A a u, 7r aIJ ) subject to the simplicity and SO(D + 1) Gaufi constraint. This way we arrive at a 
connection formulation in terms of the compact gauge group SO(D + 1) although we are con- 
sidering Lorentzian gravity. Of course, the fermionic contribution to the Hamiltonian constraint 
of Lorentzian gravity, just as in 3 + 1 dimensions, acquires correction terms as compared to its 
Euclidean counterpart which in part is due to switching from Lorentzian to Euclidean 7 matri- 
ces. Yet, these corrections are not as cumbersome as one might expect because the Lorentzian 
Clifford algebra differs from the Euclidean one just by a factor of i in front of 7 . 

After having obtained the fermionic contributions to the classical constraints we quantise 
them using standard methods [I] and using the representation [12]. Our formulation provides 
access to the quantisation techniques developed in [13^ [T31 [T5| [T6l [T7] and allows for a rigorous 
construction of a quantum theory. Ideas on how to quantise the simplicity constraint, which has 
been shown to be anomalous in [181 E] are supplied in [19] . 

2 Canonical Analysis of Lorentzian Gravity coupled to Dirac 
Fermions 

As opposed to pure gravity where, in the end, it does not matter whether one starts with a 
first or a second order formulation of the theory, this choice results in inequivalent theories 
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when dealing with fermions. The reason for this is that the torsion freeness condition which one 
derives when starting with first order general relativity, as e.g. done in our companion paper 
[2], is modified by a term quadratic in the fermions, thus resulting in a non- vanishing torsion. 
At the end of the canonical analysis, one arrives at the same set of variables, but, after solving 
the equations of motion for the torsion part of the connection, one obtains more interaction 
terms, most prominently four-fermion interactions, which are not present in the theory when 
starting with a second order formulation. To the best of the authors' knowledge, it is unclear 
which type of action should be preferred on physical grounds. The second order variant leads to 
less interaction terms and could thus be preferred by demanding simplicity. On the other hand, 
when deriving the Ashtekar-Barbero variables from the Hoist action, one deals in a first order 
framework and one could thus consider it more natural to choose this route. In this paper, we 
will choose the first order approach since the results of the canonical analysis of our companion 
paper [2J can be nicely used in order to deal with the torsion terms. For further literature on 
this topic, we refer to [10| 14"]. 

We start with the first order action 

S G + F = -J^ d D+l x (hxPe^F^M) + l -Ve<tfD^ - '-D^etf*) . (2.1) 

A4 denotes the space-time manifold of topology IR x a. Ai foliates into hypersurfaces S$ := 
X t (a), where X t : a — > M is an embedding defined by X t (a) = X(t,a). x a ,a,b,c = 1, . . . , D 
are local coordinates on a. The signature of the space-time metric is given by (— , +,..., +). 
e^ 1 denotes the vielbein and F^ u jj := d^A u jj — d u A^jj + L4 M ,A„]/j is the field strength of 
the SO(l,.D) connection A^jj. \& denotes a Grafimann valued Dirac spinor. Additionally, 
= + |A a /jS /J *, S /J = -|[7 7 ,7 J ] and * = <5 1 7°. The properties of the 7 matrices 
are summarised in appendix [AJ The gravitational part of this action has been analysed in [IJ[2], 
we will therefore concentrate on the fermionic part. 

The split in space and time is performed analogously to the D = 3 case. We split the 
time-evolution vector field into 



T M = Nn^ + N>*, = 0. (2.2) 

N is called the lapse function and the shift vector field, is the unit future pointing vector 
field normal to the spatial slices Ej, i.e. n^n^ = — 1 and n^d^Xt = 0. We use 

8»v = {g^u + n^ny) - n»n v =: q^ v - n»n u (2.3) 

to project the vielbein as 

e ^ = qi* ue » - e^r^ =: "e^ - . (2.4) 

We choose the time gauge prior to the canonical analysis (the time gauge is a canonical 
gauge, see, for instance, [20]) by setting n 1 = <5q. The action can now be rewritten ad3 

S g + f = J dt J d D x(EfK i a + i(^)(^)--N'H-N a n a -X ij G^ -(A* + (2.5) 

X A proper canonical analysis would work in terms of the momentum % aIJ conjugate to A a ij. However, 
here we use a shortcut and use the fact that the corresponding simplicity constraint [U [2] imposes ir aIJ = 
2n I[ E aJ \ n J nj = -1, mE aI = 0. 
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where 



1 1 _,„,, . „ 11 



tr.fr. 



H = -^ q R + —E^E^K ai K b] + —-K^-F^ m K m 
i 1 —^^r, i 1 



+ - — yq^EfY Va(yqV) - - — V a (^)Eff(^) 
-±(^)^(^)K ai E^ - l±^ffSE% [l\ ^} mm%f; (2.6) 

H a = -2i^ V [a i^. + ^*)t^ (2.7) 

G ij = 2K a i E a ^ - (^g*)tE«(^g*), (2.8) 
= K%E a \ (2.9) 

and small Latin indices i,j,k,... = 1, . . . , D are internal indices in the time gauge. F a ^' bkl and 
the derivation of the symplectic structure have been described in [H [2]. The Gaufi constraint has 
been split into its rotational part G %3 and its boost part G l . A™ := —T^A^ij and A.; = Ajo- All 
terms proportional to K^-E a ^ not belonging to the boost part of the Gaufi constraint have been 
written as *{G l . Further, E ai = <J(pe ai and we have decomposed A a u = T a ij + 2n^K a \j] + K a ij . 
The bar notation K a %j which was already used in our companions papers [IJ [2] means that the 
internal indices are orthogonal on n 1 . In the time gauge, this is equivalent of having only small 
latin indices running from 1 to D. The splits in trace and trace free parts are done with respect 
to the vielbein, e.g. = ^K^E bk E a i. The boost part of the Gaufi constraint does not 

acquire a fermionic part because of the cancellation S° l + (£ 0l )t = 0. The Dirac spinors in the 
above equations appear only as half-densities, i.e. ^fq^> . Since the symplectic structure tells us 
that these half-densities are the natural canonical variables, we will abuse notation and denote 
by ^ from now on a half-density. The importance of using half densities stems from the simple 
form of the symplectic structure. Otherwise, the connection would acquire a complex part [21] 
and the techniques introduced in [131 13 EH E] would not be accessible. 

In order to facilitate the canonical analysis, we will employ the equations of motion for Aj and 
K%f- at the La grangian level. Their solutions translate directly to a purely canonical treatment 
as one can check. Variation of the Lagrangian with respect to Aj sets the boost part of the Gaufi 
constraint to zero. Variation with respect to K^ r ' yields 

^• = ^*<{7 m ,^'}*, (2.10) 

which we use to eliminate K^i^' in H. 

Next, we perform the Legendre transform, yielding the constraints 

H = l^R+^-E^E^KaiK^ 

i 1 -=-„„ „•„.,. i 1 



+-— *E?f V a * - -— V a *E?7** 
2\/9 2 y/Q 

-IvWv&Vij* + ^{l k , ^} **{7fe, Etf} *, (2.H) 

H a = -2^'V [o Jf 6]i + ^t Va «r-I(V a *)t*, (2.12) 

G ij = 2K[ i E a ^ - *ts<i$ j (2.13) 

where V a is the covariant derivative associated to the spin connection T a ij^ &s well as the non- 
vanishing (generalised) Poisson (anti-) brackets [20] 

{E ai (x),K bj (y)} = 8 D (x-y)5^i and {* a {x), -M\{y)} = -5 D {x - y)5 a p . (2.14) 
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A term proportional to the Gaufi constraint has been omitted in H and T-L a - 
We define the generator of spatial diffeomorphisms 



which acts as 



1 



T aij G^ = -E b ^d a K bj + d b (E b iK aj ) + ±&d a * - ^d a tf)V, 



[E a \n b [N b ]] 

K au Hb[N b }} 
{y,H b [N b ]} 
[^,n b [N b ]} 



N b d b E ai + (d b N b )E ai - (d b N a )E H 
N b d b K ai + (d a N b )K bu 
N b d b V+^(d a N a )V, 



by Lie derivatives. The GauB constraint acts as 

1 



E 



1 



1.1 J 



\ % jE a \ 
X^K 



*r{ 7 fe ,[A,S]^} + {A fem7 TO ,E^}) 



(2.15) 



(L> a *)t,iG^'[A 
{*'>:''•''*. A,/ 

*{7 fe ,S^}*,i^[A iy ] 

We therefore conclude that the algebra of the diffeomorphism and GauB constraints closes 
and that they both Poisson-commute with the Hamiltonian constraint, at least weakly. 

Thus we are left with checking the Poisson bracket of two Hamiltonian constraints. We 
split H = Hgrav + %2F + %4F into the purely gravitational part, a part containing two fermions 
and a part containing the four-fermion terms and define V a := Md a N — Nd a M as well as 
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V a b := (d a M)(dbN) — (d b M){d a N). The non-vanishing Poisson brackets are given as 

{n gTav [M],n grav [N]} = J d D xlv a q ab (-2E^V [b K c]j ) (2.29) 



E ai E bj ^ \ 

Vab K a[i E j] , 



{%F[M],% P [iV]} = Jj D x (v a q ab ^W a y- l -(V a ^^ 




{H 2 f[M],H4f[N]} + {H&[M\,Hzf[N]} = j dPx f-l-V a ^{E^ k ^jW^.y 



and sum up to 

{H[M],H[N}} = jj D x (v a q ab H b + Vab^-G'A . (2.30) 



The constraints are therefore consistent and the canonical analysis ends here. 

3 Phase Space Extension 

In [TJ [2] it was shown that the extension of the ADM phase space (q a b,P ab ) to the extended 
phase space (A a jj, ir aIJ ) subject to Gaufi and simplicity constraint is equivalent to the ADM 
phase space. Moreover, this is possible using SO(D + 1) as the structure group while considering 
Lorentzian gravity. Since spinors can only be coupled to vielbeins, we have to construct a 
transformation from (E at ,K a i) to (A a jj, 7r aIJ ). The calculation turns out to be very similar to 
the one described in [H [2], we therefore only give the result and comment on some peculiarities. 
The explicit construction is given by 

E aI = Cfj I JK aKJ n K , K aI = (f] I J (A-T) aK jn K , (3.1) 

where ff J = rj IJ — Qn 1 ^ ~ rf J — {n aK1 ^aK J — Cv IJ )i ^aU is the hybrid connection of 
E aI (see [TJ [2] for details) and r/ 00 = r\ % i = 8^. The peculiarity of these expressions is the 
appearance of n 1 , which can only be directly (that is, without non-polynomial terms except for 
y/q) expressed in terms of it aIJ for D + 1 odd. For general D, we only have access to n Ir n J and 
then can define ±n^ through ±n J = V n / n / sgn(n°n^) (no summation understood here and one 
substitutes for n^n 1 under the square root the expression for n 1 ^ at I = J). Fortunately, we 
can avoid to make use of this explicit square root expression by invoking the following trick: 
Ultimately the non- vanishing Poisson bracket involving n is of the form n {A a jj,n K }. Since 
n K riK ~ C modulo simplicity constraint we have n J {A a jj,n K }riK ~ 0. To see this, notice that 
the simplicity constraint reads Sjjj- = ^ijklm^ 1 ' 7 7T (see [H[2] for details). It follows 

n J {A aIJ ,S£} = 2n J e IJKLl j5^ KL « 

on the constraint surface TT aIJ = 2n^E aJ \ It follows n J {A a u,n K } « n J {A a u ,n L }f)f = 
—n J {A a jj,f]^}n L . However, {A a jj, fj^} = —^{A a jj,n K ni,} and n K riL can be expressed un- 
ambiguously as above in terms of 7r aIJ . In order to compute the brackets between E aI ,K a j 
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one then just hast to carefully insert the definition of njnj in terms of 7r a/J ' . The only term 
which cannot easily be seen to vanish by algebraic manipulations alone occurs in the bracket 
{K a j, K a j} and is of the form 

f)fn L f)f n N {[A - T] aIJ , [A - T) hKL } = -fjf n L fjfn N [{[A aIJ) Y hKL } - {A bKL ,T aIJ }]. 

This term vanishes due to the weak integrability (modulo simplicity constraint) of the hybrid 
connection T a jj and by using the trick mentioned above. See pQ for more details. 

After a tedious calculation, the Poisson brackets of E aI and K a j expressed as functions of 
A a jj and 7r aIJ are given by 

{E aI (x),E bJ (y)} = 0, {K aI (x),K b j(y)}*0, {E aI (x), K bJ {y)} « -(5 D (x - y)5f5 I J . (3.2) 

modulo simplicity constraint. 

The only task left to do is to write down a Hamiltonian theory in the variables A a jj and 
Tr bKL with internal gauge group SO(D + 1) which reduces to the theory derived in the previous 

section on the constraint surface S£ = njG IJ = 0. The basic idea is to first derive a Hamilto- 

M 1 

nian formulation of Euclidean gravity coupled to fermions and then to adjust the Hamiltonian 
constraint to mimic Lorentzian gravity. The reason why the procedure introduced in [TJ [2] gen- 
eralises nicely to Dirac fermions is the strong resemblance of the Clifford algebras, which differ 
only by factors of i for different signatures and the Euclidean signature of the internal gauge 
group which ensures that E is a Hermitian matrix the Euclidean case. This requires care at 
several places, e.g. the cancellation S 0i + (X°*)t = from the boost part of the Lorentzian Gaufi 
constraint is no longer present. In order to derive the Euclidean constraints, we start as in the 
previous section with the action (|2.ip and perform a D + 1 decomposition. We replace 7 with 
n/7 , which reduces to 7 in the time gauge. We note that the object is not a Lorentz 
scalar any more when using Euclidean signature, because the 7 inherent in tyty is needed in 
order to maintain invariance under boosts which are generated by the anti-Hermitian £°\ In 
Euclidean signature the boost generator is also Hermitian and thus ty'ty rather than is now 
the appropriate Euclidean scalar to be used while vpty^vl/ [ s a Euclidean covariant vector with 
index /. The substitution 7 — > nij 1 is therefore natural for Euclidean signature and allows 
for the construction of a manifestly SO(D + 1) gauge invariant theory. We use the additional 
n 1 in the action to form Tr' aIJ = 2n^E a ^ (see [U [2] for details) and introduce the simplicity 
constraint in order to replace ir' aIJ by 7r aIJ . The Euclidean Hamiltonian theory is then given 
by the constraints 



(3.3) 



H E = \ aIK TT bJ K F abI j+{U^ aIJ ^jjD a ^ + CC^ 

U E a = l -K bIJ F abI j+ % -^D a ^> - % -{D a V?V, (3.4) 

G 1 ^ = D a 7r aIJ - ^Z IJ ^, (3.5) 

oab ^ , „alj bK L fn c\ 

b TT = ~a IJKLM 71 > ( 3 - 6 ) 



M — 4 C IJKLM 

and the Poisson bracket 



{A aIJ (x), 7r bKL (y)} = 5 D (x- y)5 b (5 I K S J L - S 1 ^), = -S D (x - y)8$. (3.7) 

The task of "Lorentzifying" the gravitational part of H E has been addressed in [HH]. For the 
fermionic part, we observe that we should add a factor of i in front of the fermionic term in 
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order to compensate for 7^ = ij^ and denote the changed constraint by 7~t F y The Hamiltonian 
constraint now reduces to 

-^L-tf*E ai K ai - -L^l^^^ + 9a ^«rt s «^ + 0(i^-).( 3 .8) 
The terms proportional to K^j' can be dealt with using ideas from gauge unfixing. We calculate 

:= fl£[JVcg] + I^S (3-9) 

and see that the gravitational constraint DjL = —e IJKL j^K cIJ (n^ KN D c ir b ' L n) ~ now re- 
ceives a fermionic contribution pD. F4%%£, the Dirac matrix introduced in [2], however remains 
unchanged since pD Poisson-commutes with the simplicity constraint and gauge unfixing works 
as before. Next to compensating the terms proportional to i^ij > gauge unfixing also produces 
a four-fermion term, which we have to subtract again in order to build the correct Lorentzian 
Hamiltonian constraint. 

Comparison with the previous section leads to the following correction terms: 

n L = nf l) + ^E^E^K aI K u -l G D^(F^)Y b f dG D§ 
~2 m ' ' ab cd F n ~ 2 F M ' ' ab cd G w 



2^/g ^2^/gD-l 01 

-a a ^^ f Sl*J + ^=^£ a/ ^ a/ + ^{7*, E^}^{ 7fe , Ey}*. (3.10) 

This Hamiltonian has to be rewritten in terms of A a jj and n bKL only, desirably as simple 
as possible regarding the quantisation. We propose the Hamiltonian 

n L = ]^ aIK 7: b \F aM j+^U>\ aU T 1I jD a ^ + CCj 

+l-E^E b ^K aI K bJ - 1 - G D^ (F-^JSonfi 

n ab (t?-1\WN ncd 1 nab (j?-l\TT~N ncd 
~2 G M ' > ab cd F N ~ 2 F A7 {* ) ab cd G^Tv 

^E IJ ^Eu^f + — L 3 ^ ~ A ^T l IK ^T lJK ^n I n J (3.11) 



2^q 2^qD-l 
-d a (^-tfE IJ *) + ^-^E ai K ai - ^h [I l J l K l L] ^n L n M ^ 1[nJlKlM] ^ 

for quantisation, although we are well aware of the fact that other choices might lead to equally 
reasonable classical starting points. We shift the problem of choosing the "correct" Hamiltonian 
constraint to the semiclassical analysis. The expressions for rJnj and E aI K bI were derived in 
[TJ [2] and all 7-matrices appearing are those for Euclidean signature. 
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4 Kinematical Hilbert Space for Fermions 



The construction of the kinematical Hilbert space for fermions was discussed in |12| . Results 
obtained there apply for the case at hand, so we only give a short summary. It is crucial to work 
with half-densitised fermionic variables ^ for what follows, as was stressed in \12\ . 

Faithful implementation of the reality conditions enforces the use of a representation in which 
the objects 



t (x) 



/ d D y^6&yj ■= lim / d D y*^pfi* a (y) (4.1) 



become densely defined multiplication operators. Their adjoints 9 become derivative operators. 
Here, a = l,...,n := 2^ D+1 ^ 2 ^ (|_.J denotes the integer part of .) and Xe(x,y) denotes the 
characteristic function of a box of Lebesgue measure e D centered at x. In the above equation, 
the half-densities ^ are "dedensitised" using the (^-distribution, which is a scalar in one of 
its arguments and a density in the other. Thus, the variables 9 are Grafimann-valued scalar 
quantities, which is important for diffeomorphism invariance [12]. In calculations it is understood 
that the e — > limit is performed after the manipulation under consideration is performed. 

The variables a (y) coordinatise together with their conjugates the superspace S v at the 
point v. The quantum configuration space is the uncountable direct product S := Y\ vea S v . In 
order to define an inner product on S, it turns out to be sufficient to define an inner product 
on each S v coming from a probability measure. The "measure" on S v is a modified form of the 
Berezin symbolic integral [22] 



dm (9, 9) = dOdde 60 and dm v = <g>" =1 dm (9 a (v),6 a (v)) , (4.2) 

which has the additional property of being positive on holomorphic functions (those which only 
depend on the 9 a and not on 9 a ). Since the 9 are Grafimann variables and thus anti-commute, 
any product of more than n of these variables will vanish. The vector space of monomials of 
order k is n\/k\{n — fc)!-dimensional (0 < k < n) and the full vector space Q v built from all 
monomials has dimension 2™. The full fermionic Hilbert space is a space of holomorphic square 
integrable functions on S with respect to djjLp 

% F = L 2 (S,dfi F ) = (g) 1 , eo -L 2 (S v ,dm v ) . (4.3) 

When restricted to a point v, the inner product can be seen to coincide with the standard 
inner product on Q v when viewed as a vector space of exterior forms of maximal degree D + 1. 
For a more complete treatment, the reader is referred to [12] where it is shown that the fermion 
measure dfiF is gauge and diffeomorphism invariant and that the reality conditions 9 a = —iit a 
are faithfully implemented in the inner product. 



5 Implementation of the Hamiltonian Constraint Operator 

The quantisation of the purely gravitational Hamiltonian constraint in dimensions D + 1 > 3 
was described in [3] and we will not repeat it here. The quantisation of fermionic degrees of 
freedom was described in detail in [U [12], which we assume the reader to be familiar with. 
Next to an explicit example, we will only provide a toolkit to quantise the fermionic part of the 
Hamiltonian constraint operator as writing down the explicit terms is rather laborious. 

Quantisation of the 9 variables is performed by promoting 9 a to a multiplication operator 

and (9^)P = — h-§Q-, where L indicates the left derivative. The explicit quantisation follows the 
(extended) toolkit of [3]: 
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(1) Choose a triangulation T(7, e) of the spatial slice a adapted to the graph 7. 

(2) Use classical identities in order to express the Hamiltonian constraint in terms of connec- 
tions A a jj, volumes V(x, e) and Euclidean Hamiltonian constraints "He(x,€). 

(3) Replace all phase space variables by their corresponding regulated quantities. 

(4) Instead of the the integration d D x, put a sum ^^^^ over all the vertices v of the 
graph 7. 

(5) For every spatial e-symbol, put a sum -J^y Y2v(A)=v over an -D-simplices having v as 
a vertex. The holonomies associated with the e-symbol are evaluated along the edges 
spanning A. 

(6) Substitute the generalised Poisson (anti-)brackets by i times the (anti-)commutator of the 
corresponding operators, i.e. the multiplication operator h e , the volume operator V, the 
multiplication operator 9 a and the derivation operator —h^-. 

The kinetic fermionic part of the Hamiltonian constraint operator is a bit more involved 
since it contains a derivative. Following [12], we explicitly get 

ffDta, ktoW/, (5.1) 

= (f E gig) E '(^(^(^((WMAKe))-*))^^.^,, 

where by e (. . .) we mean the regulated quantity with the Poisson brackets substituted by i/h 
times the commutator of the corresponding operators. The Hermitian conjugation operation 
H.C. is meant with respect to the inner product on the Hilbert space. 

Due to its length, we refrain from writing down the complete Hamiltonian constraint operator 
which can be easily done when following the quantisation recipe. We remark that we could split 
the Dirac fermions for D + 1 even into left- and right-handed parts, however, the presentation 
does not benefit from this. Details are supplied in [23]. The quantisation ambiguities from loop 
quantum gravity are also present when considering fermions and, as usual, we shift this problem 
to the semiclassical limit. 

6 Outline of the Quantum Simplicity Constraints 

In this section, we briefly outline the implementation of the quantum simplicity constraint. As 
shown in [3], the quadratic simplicity constraint can be implemented on edges anomaly free 
and the result coincides with the spin foam treatment by Freidel, Krasnov and Puzio [24j . The 
vertex constraints are second class when using holonomies and fluxes as basic variables and thus 
require a more elaborate treatment. One possibility is the master constraint approach discussed 
in [3]. Another approach is discussed in [19| . which only quantises a maximally commuting 
subset of vertex simplicity constraints and thus closely resembles the gauge unfixing procedure. 
A benchmark which a successful proposal for the simplicity constraints in the canonical theory 
have to satisfy is anomaly- freeness with respect to the action of the Hamiltonian constraint. 
More precisely, the Hamiltonian constraint should leave the kernel of the simplicity constraints 
invariant. At present, none of the proposed treatments is entirely satisfactory and they deserve 
further research. We refer to our companion papers (3j [19] for details. 
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7 Conclusion 



In this paper we extended the new connection formulation of General Relativity in all spacetime 
dimensions D + l>3[TJ[21[3]to arbitrary gauge boson (for compact structure groups) and stan- 
dard fermionic couplings. We did not discuss Higgs fields or other scalars but remark that the 
strategy for background independent quantisation of Higgs field coupled to General Relativity 
developed for 3 + 1 dimensions [1] straightforwardly generalises to higher dimensions. Fermionic 
fields deserve special attention when we switch from the Lorentzian to the Euclidean Clifford 
algebra which is a necessary step if we want to maintain a compact structure group for Gen- 
eral Relativity in order to have access to background independent Hilbert space representations 
|17j . In this paper we showed that while there arise corresponding correction terms, there is no 
problem in principle in performing this switch and one can apply the quantisation techniques 
developed for 3 + 1 dimensions. With this work, the stage is set to study Supergravity theories, 
which we pursue in [231 EH] • 
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A Gamma Matrices 

The properties of the gamma matrices can be found in most textbooks on quantum field theory, 
see, for instance, [26J. Their basic property is the Clifford algebra 

{ 7 / ,7 J } = 2 ?7 /J , (A.l) 

where rj IJ is the flat Minkowski metric of a spacetime with signature (p,q). From this relation 
alone, one deduces, 

[ E IJ ys-] = -rfrfK + i^JJK (A.2) 

and 

i [S /J , S^ L ] = n LJ E KI - V LI X KJ + V JK Z IL - V IK Z JL , (A.3) 

where T, IJ := — \ [7 7 ,7 J ]- X /J thus constitutes a representation of the Lie algebra so(p, q) on 
spinor space. 

Furthermore, the expression {7^, T, IJ } = —i*y[ K *y J 7^ is completely antisymmetric in /, J, K. 

It is noteworthy that S /J is a Hermitian matrix for Euclidean signature. In general, = 
r] n r] JJ T, IJ , which becomes important when dealing with Lorentzian signature, i.e. the boost 
part of the Gaufi constraint is purely rotational as S°* + = 0. 

Explicit representations of the gamma matrices exist for all dimensions D + 1 > 2, see, for 
instance, |27j , or [28J. A generalisation of left- and right-handed spinors exists for D + 1 even 
and is spelled out in our companion paper [23]. 
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